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SECTION  1 


INTRODUCTION 


This  volume  contains  instructions  which  enable  a  non-programmer 
to  use  the  Spiral  Decay  Program  (SPIRDECA)  under  the  control  of  the  B-3 
executive  on  the  Philco  212  computer.  The  program  has  been  designed  for  use 
in  the  SPACETRACK  center  of  NORAD . 

The  theoretical  details  have  been  presented  in  Volume  I  of  the 
report.  Programming  details  will  be  presented  in  Volume  III. 

The  diverse  uses  and  many  options  of  the  program  demand  some 
attention  on  the  part  of  the  analyst  in  submitting  cases.  It  is  suggested 
that  the  input  deck  be  reviewed  immediately  before  submission  with  special 
attention  to  the  parameter  cards.  For  the  utmost  precision,  it  is  necessary 
to  use  the  best  values  of  solar  activity  indices  (supplied  by  the  Air  Weather 
Service)  and  geopotential  coefficients  (P-cards  6,  7  and  8).  A  few  seconds 
forethought  may  save  much  computer  and  analyst  time.  If  it  is  necessary  to 
submit  a  series  of  runs  with  very  similar  decks,  use  should  be  made  of  the 
Remark  cards  permitted  by  the  Executive  System.  These  cards  have  the  letters 
REM  in  columns  17-19,  columns  20-24  blank  and  the  alphanumeric  remarks 
following.  Remark  cards  may  be  placed  immediately  before  the  SPS  JOB  card. 

During  the  training  and  experimentation  at  SPACETRACK,  analysts 
have  found  operation  of  the  program  to  be  comparatively  easy.  Many  of  the 
operational  features  of  the  program  were  suggested  by  personnel  of  the  1st 
Aerospace  Surveillance  Squadron. 
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SECTION  2 


SPIRAL  DECAY  PROGRAM  OPERATING  INSTRUCTIONS 


The  B-3  version  of  the  Spiral  Decay  Program  (SPIRDECA)  has 
been  modified  to  eliminate  OBSWGT  and  the  OBSWGT  control  card  (P  Card  0) , 
and  to  accept  984,  rather  than  400,  observations.  The  program  may  be 
used  in  three  modes:  (1)  Differential  Correction,  (2)  Prediction  or 
(3)  a  combined  Differential  Correction  and  Prediction. 

2,1  INPUT/OUTPUT  INSTRUCTIONS 

Examples  of  the  job  deck  setups  for  each  mode  of  operation 
are  shown  in  Figures  1-3. 
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Logical  tapes  0,  7,  12  are  optional  and  should  be  used  for  the 
following: 


Logical  0  -  if  the  number  of  observations  exceeds 

492 

Logical  7  -  to  override  or  add  to  the  assembled 

weights  and  biases 

Logical  12  -  if  a  binary  ephemeris  tape  is  requested. 

The  program  output  and  a  dump,  if  taken,  is  printed  from 
logical  tape  11  on  Data  Select  1.  Cards  are  punched  on  Data  Select  2. 

The  binary  ephemeris  tape  (Logical  12)  may  be  used  as  input 
to  two  programs,  XYZLA  and  PUNCH,  to  generate  look  angles. 

SPIRDECA  may  also  be  interfaced  with  the  GIPAR  program.  The 
input  decks  must  be  stacked  and  run  together.  Five  points  may  be  passed 
to  GIPAR  by  use  of  P  Card  10  (Figure  18  ) . 
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FIGURE  1 .  SPIRDECA  INPUT  DECK  FOR  DIFFERENTIAL  CORRECTION  RUN 
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FIGURE  2 .  SPIRDECA  INPUT  DECK  FOR  PREDICTION  RUN 
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FIGURE  3  .  SPIRDECA  INPUT  DECK  FOR  COMBINED  DC  AND  PREDICTION 
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FIGURE  8.  END  SCHEDULE  TAPE  CARD 
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2.2  PARAMETER  CARDS 


Card  formats  for  the  input  parameter  cards  are  included 
in  Figures  9-18. 

a .  P  Cards  1  and  2 


Elements  may  be  input  in  two  formats  (r,  r  or  L,  a,  h)  on 
P  cards  1  and  2.  When  this  is  done,  a  6  card  element  set  is  not 
required  as  input. 

b.  P  Card  3 

The  information  necessary  to  control  the  differential 
correction  is  on  P  card  3 

c  .  P  Card  4 


The  control  information  to  run  a  prediction  is  on  P  card  4. 

d.  P  Card  5 


This  is  on  optional  card  which  contains  error  criteria  for 
the  Adams-Bashforth  integration.  It  should  be  used  only  to  override 
the  assembled  values,  which  are  currently  10”8 ,  for  each  element. 

e.  P  Card  6 


P  card  6  is  an  optional  card  and  contains  the  zonal 
harmonic  terms  to  be  used  in  the  bulge  perturbations.  It  is  necessary 
only  when  one  wants. to  change  the  assembled  values,  which  are: 

J2  =  1082.50  x  10-6 

J3  =  -2.553  x  10-6 

J.  =  -1.807  x  10-6 

4 

J5  =  -8.3  x  10-8 

f .  P  Cards  7  and  8 


If  tesseral  harmonics  are  needed  in  the  bulge  perturbations, 
the  coefficients  (S  and  C  )  must  be  entered  on  P  cards  7  and  8. 

Both  cards  must  be  ?n  the  input  deck,  if  this  option  is  used,  even  if 
all  the  fields  are  not  filled  in. 
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g.  P  Card  9 


This  optional  card  controls  the  printed  output  of  the 
prediction  ephemeris.  Only  the  intervals  specified  will  be  printed. 

h.  P  Card  10 


This  optional  control  card  is  necessary  only  when  GIPAR  will 
be  run  immediately  after  SPIRDECA.  The  times  specified  will  be  given 
to  GIPAR  and  also  printed  if  the  print  option  is  set. 
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FIGURE  9.  ELEMENT  CARD  1 
(P  CARD  1) 
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FIGURE  10.  ELEMENT  CARD  2 
(P  CARD  2) 
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14 
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1  Print  input  obs.  (0  =  No,  1  =  Yes) 

2  Weight  flag  (0  =  Use  assembled  weights,  1  =  Read  from  tape) 

3  Print  weights  and  biases  (0  =  No,  1  =  Yes) 

4-8  Not  Used 

9  New  Epoch  mode  (0  =  by  Time,  1  =  by  Rev.,  2  =  Time  of  last  obs.) 

10-24  New  Epoch  time  (YYDDDHHMMSS .SSS)  or  Rev.  No.  (NNNNN) 

25  Not  Used 

26-32  Elements  to  correct  (n,  a  ,  a  ,  U  ,  f),  i,  B) 

33  Max.  number  of  differential  corrections 

34  Correct  n  only  on  1st  pass  (0  =  No,  1  =  Yes) 

35  Aq  check  (0  =  No,  1  =  Yes) 

36-39  Maximum  Aq  (km) 

40-42  1st  pass  rejection  (km)  for  range,  azimuth,  elevation 

43-45  1st  pass  rejection  (km/sec)  for  range  rate 

46-48  RMS  multiplier  (Recommended  value:  4) 


FIGURE  11.  DIFFERENTIAL  CORRECTION  CONTROL  CARD 

(P  CARD  3) 
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FIGURE  11.  DIFFERENTIAL  CORRECTION  CONTROL  CARD  (continued) 

(P  CARD  3) 
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25  Not  Used 

26  Not  Used 

27  Prediction  output  flag  (0  =  none,  1  =  t,  £,  r; 

2  =  t,  0,  X,  h;  3  =  both) 

28-40  Not  Used 

41  Binary  tape  output  flag  (Ephemeris  tape)  0  =  No,  1  =  Yes) 

42-45  At  (min)  for  Ephemeris  tape 

46-78  Not  Used 

79  Card  Number  =  4 

80  Card  Type  =  P 


FIGURE  12.  PREDICTION  CONTROL  CARD 

(P  CARD  4) 
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FIGURE  14.  ZONAL  COEFFICIENT  CARD 

(P  CARD  6) 
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FIGURE  16.  TESSERAL  COEFFICIENT  CARD  (2  of  2) 

(P  CARD  8) 
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FIGURE  18.  GIPAR  CONTROL  CARD 
(P  CARD  10) 


24 


2.3  AUXILIARY  WEIGHT  TAPE 


To  change  or  add  station  weights  or  biases  in  the  SPIRDECA 
program,  an  auxiliary  input  tape  must  be  mounted  on  logical  tape 
unit  7.  The  deck  setup  for  this  tape  is  shown  in  Figure  19.  The  first 
card  is  an  identifying  header  and  the  last  card  an  end  card  (see 
Figure  20) .  Between  these  cards  are  the  bias-weight  input  cards 
(Figure  21)  . 

This  input  tape  is  required  only  when  new  weights  or  biases 
are  desired.  Storage  in  the  program  is  limited  to  30  sensors. 
Therefore,  the  number  of  distinct  sensors  appearing  in  the  assembled 
program  plus  those  on  the  auxiliary  weight  tape  must  not  exceed  30. 

The  data  on  the  cards  override  those  in  the  assembled  program.  When 
the  auxiliary  input  tape  is  to  be  used,  column  2  on  P  card  3  must 
have  a  one-punch  (1)  to  inform  the  program  of  the  existence  of  the 
weight  tape. 
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(TO  BE  USED  AS  INPUT  ON  LOGICAL  TAPE  7) 


/  ENDSIGKA 

/  (Weighting  Cards) 

/  70WEIGHT  ^ 

FIGURE  19.  INPUT  DECK  SETUP  TO  USE 

WEIGHTS  OR  BIASES  IN  OBSWGT 
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FIGURE  21.  BIAS  -  WEIGHT  CARD  (EDITOR  AND  B-3  MASTER) 
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SECTION  3 


CALIB  OPERATING  INSTRUCTIONS 


The  B-3  version  of  the  Sensor  Calibration  Program  (CALIBA)  has  been 
modified  to  eliminate  OBSWGTA  and  the  OBSWGTA  control  card;  and  to  accept  984, 
rather  than  400,  observations, 

3.1  INPUT/OUTPUT  INSTRUCTIONS 

a  .  Input  Deck  Setup 

The  CALIBA  program  requires  four  parameter  cards  in  the  standard 
SCHED  tape  job  deck. 

P  Cards  1  and  2  introduce  input  elements  and  epoch,  and  are  identical 
to  those  used  in  SPIRDECA  (Section  2,  Figures  9  and  10).  P  Card  0  format  is 
given  in  Figure  24;  P  Card  3  format  is  given  in  Figure  25. 

The  P  Card  images  are  automatically  printed  by  the  SPIRDECA  program 
to  simplify  the  interface  between  these  programs. 
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b .  Tape  Setup 


Logical  0 

Scratch 

WR 

1 

B-3  Master 

2 

Input 

WR 

4 

SEAI 

WR 

11 

Output  Tape 

WR 

Logical  tape  0  is  optional  and  should  be  used  if  the  number  of 
observations  exceeds  492. 

c .  Output 

The  program  output  tape  is  printed  in  Data  Select  1  for  program 
results  and  dumps,  if  taken. 

The  CALIBA  program  has  no  teletype  or  punched  card  output. 
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NOTE:  All  card  formats  referenced  in  this  figure  may  be  found  in  Section  2 

with  the  exception  of  the  SPSJOB  Card  (Figure  23),  P  Card  0  (Figure  24), 
and  P  Card  3  (Figure  25) . 
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FIGURE  22.  CALIBA  INPUT  DECK  SETUP 
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FIGURE  23.  SPS JOB  CARD 
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3.2  PARAMETER  CARDS 


The  card  formats  for  the  input  parameter  cards  included  in  this 
section  are: 

a .  P  Card  0 

Sensor  weight  information. 

b.  P  Card  3 

Information  necessary  to  control  the  sensor  calibration. 
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3.3 


OBSERVATION  DATA 


The  observation  data  cards  used  in  the  CALIBA  program  must 
(1)  be  from  the  same  sensor ,  (2)  be  from  the  same  satellite,  and  (3)  be 
of  the  same  "observation  type" . 
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SECTION  4 


SUGGESTED  OPERATING  PROCEDURES 
SPIRAL  DECAY  PROGRAM 


These  guidelines  for  the  application  of  the  Spiral  Decay  Program 
are  based  upon  several  months  of  program  application  in  the  ADC  and  496L 
computation  environments.  Since  an  improved  SOP  will  only  emerge  from 
continued  application  and  documentation  of  experience,  suggestions  for 
improvement  in  procedures,  documentation  or  program  output  should  be 
directed  to  Aeronutronic  personnel  at  Colorado  Springs  and  at  Bedford. 

The  majority  of  the  guidelines  provided  herein  apply  to  Spiral 
Decay  exercises  with  ADC  data.  Some  generalizations  to  higher  altitude 
satellites  and  to  SCF  data  (PRELORT/VERLORT)  are  inferred. 

The  most  frequent  source  of  difficulty  in  using  the  Spiral  Decay 
Program  involves  input  card  errors.  A  high  degree  of  flexibility  has  been 
built  into  the  program  and  control  of  available  program  options  necessitates 
fairly  complicated  parameter  card  formats.  The  few  minutes  required  of  the 
analyst  to  recheck  the  input  deck  may  avoid  many  false  starts  on  the  computer. 

4.1  INITIALIZATION 

Initialization  may  be  accomplished  from  P  cards  l^and  2,  with 
either  L,  a,  h  elements  or  position  and  velocity  (x,  y,  z,  x,  y,  z)  ,  or 
with  standard  6-card  element  sets.  Where  decay  is  imminent,  say  within  10 
days,  initialization  from  6-card  elements  should  be  done  in  two  steps,  as 
follows : 


(1)  With  approximately  24  hours  of  observations  and 

with  a  6-card  element  epoch  falling  within  this  24- 
hour  period,  correct  n,  UQ  and  B  only.  For  initial 
estimates  of  B,  use  the  following: 
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Payloads  B  =  0.01  m2/kgm 

Rocket  Bodies  B  =  0.02  m2/kgm 

(2)  With  the  final  elements  provided  by  SPIRDECA  in 
the  form  of  P  cards  1  and  2,  rerun  the  same  data 
package  as  a  7 -element  fit.  (Be  sure  to  change 
the  input  option  on  SPSJOB  card) . 

An  initial  three-element  fit  is  necessary,  since  the  estimates  on  B  may  be 
grossly  in  error.  These  estimates  on  B,  incidentally,  are  characteristic 
of  tumbling  bodies,  for  stable  objects,  lower  values  are  common.  B’s  ob¬ 
served  in  experimentation  range  from  0.0034  for  object  890,  a  stable  booster- 
payload  combination,  to  0.035  for  object  707,  a  tumbling  Soviet  rocket 
booster  . 

4.2  CHOICE  OF  EPOCH 

Two  or  more  epochs  must  be  defined  in  the  course  of  setting  up 
a  Spiral  Decay  deck.  These  are  the  element  epoch  (P-card  2),  new  epoch 
(P-card  3) ,  and  final  prediction  epoch  (P-card  4) . 

Normally,  an  element  epoch  within  the  data  interval  is  preferable 
to  one  outside,  both  to  aid  convergence  and  to  avoid  unnecessary  computation 
time.  As  the  exercise  proceeds,  the  element  epoch  must  be  updated,  usually 
with  each  new  addition  of  data.  The  most  convenient  way  of  automatically 
advancing  epoch  is  to  update  the  epoch,  on  each  fit,  to  the  latest  obser¬ 
vation  time  in  the  data  package,  a  program  option  available  to  the  analyst. 

In  this  manner,  the  analyst  will  receive  updated  element  cards  for  both  the 
current  element  epoch  as  well  as  the  time  of  the  last  observation. 

The  final  prediction  epoch  controls  the  prediction  program.  If  an 
ephemeris  to  decay  is  desired,  the  final  prediction  epoch  may  be  any  time 
beyond  the  decay  time.  If,  on  the  other  hand,  only  acquisition  messages 
are  to  be  generated,  through  the  binary  ephemeris  tape  option  and  the  XYZLA 
program,  the  final  prediction  time  will  define  the  final  time  for  such 
messages. 

4.3  SELECTION  OF  DATA 

For  objects  approaching  decay,  data  distribution  has  an  important 
bearing  on  convergence.  Data  should  be  well  distributed  in  the  argument  of 
latitude,  U,  as  well  as  in  time.  (For  sensors  available  to  SPACETRACK,  U 
is  limited  to  the  northern  hemisphere,  i.e.,  0°  <  U  180°).  The  majority 
of  the  data  should  be  derived  from  trackers,  with  preference  to  those  given 
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high  weight  (e.g.,  336,  337  and  345), 


If  the  data  are  concentrated  at  the  ends  of  the  time  interval, 
with  only  scattered  data  between,  a  seven  element  fit  may  fail  to  converge. 
This  is  due  to  the  fact  that  the  program  cannot  independently  determine 
period  and  rate -of -change  of  period  from  data  limited  to  the  boundaries  of 
the  time  interval. 

The  required  data  interval  to  successfully  correct  B  depends  upon 
perigee  altitude  and  on  B.  For  the  last  day  of  lifetime,  intervals  as  short 
as  a  half-day  may  be  successful,  since  the  perigee  will  be  ^  150  km.  On 
the  other  hand,  a  perigee  altitude  of  400  km  may  require  a  two-day  span;  at 
higher  perigee  altitudes,  drag  correction  may  be  impractical  unless  (1)  B 
is  large,  i.  e.  a  balloon,  or  (2)  the  data  are  exceedingly  good,  with  long 
tracks  of  both  range  and  range -rate  data.  Since  each  case  is  highly  indiv¬ 
idual,  depending  on  the  data  span  and  distribution,  perigee  altitude,  B  and 
data  quality,  a  narrow  guideline  cannot  be  established. 

Once  the  value  of  B  has  stabilized,  6-element  fits  with  short 
( ' 12  hour)  data  spans  may  be  used.  This  technique  would  be  employed,  for 
example,  where  data  near  the  decay  revolution  were  being  processed,  since 
a  7 -element  fit  under  these  circumstances  may,  in  the  iterative  correction 
procedure,  cause  the  satellite  to  "decay"  prematurely.  No  recovery  is 
provided  in  the  program  for  this  eventuality. 

The  versions  of  SPIRDECA  on  the  B-2  and  B-3  masters  are  limited 
to  two  groups  of  492  observations  (about  2.75  inches  of  cards).  An  input 
option  must  be  changed,  depending  upon  whether  the  number  of  observations 
exceeds  492.  The  program  will  sort  these  observations,  so  their  time  order 
in  the  input  deck  is  immaterial. 

4.4  INTERPRETATION  OF  RESIDUALS 

Residuals  may  be  viewed  as  a  group,  through  the  weighted  rms ,  or 
individually. 

The  weighted  rms  reflects  the  quality  of  the  fit  against  norms 
established  for  the  sensors.  In  theory,  if  the  sensors  behave  in  the  manner 
described  statistically  with  weights  and  biases,  and  if  the  systematic  errors 
in  the  residuals  are  small  due  to  a  good  fit  to  the  orbit,  the  weighted  rms 
should  approach  unity.  Larger  values  are  more  common,  due  to  these  factors: 

(1)  Sensor  performance  is  poorer  than  that  described 
with  the  weights  and  biases. 
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(2)  Large  uncorrected  signal  components  exist  in  the 
residuals,  which  cannot  be  reduced  by  the  program. 

The  latter  situation  is  common,  for  example,  where  the  satellite  ceases 
tumbling  during  the  data  span  (i.e.  the  model  assumes  a  single  drag  para¬ 
meter  for  the  entire  data  span) ,  or  where  the  data  distribution  forces  a 
limitation  in  the  program* s  degrees  of  freedom  (i.e.  one  must  make  a  3, 

4,  5,  or  6-element  fit,  rather  than  a  7-element  fit). 

Normally  a  weighted  rms  greater  than  5  indicates  excessive  signal 
component  in  the  residuals  and  the  causes  should  be  investigated  further. 

If,  for  example,  SOI  data  indicate  that  the  tumble  rate  has  been  decreasing, 
one  may  suspect  that  stabilization  has  taken  place,  particularly  if  the 
ballistic  coefficient  B  has  decreased  more  than  5  percent  from  the  pre¬ 
vious  fit.  In  such  a  case,  one  would  delete  more  old  data,  making  the 
next  7-element  fit  with,  say,  an  18-hour  data  interval. 

One  of  the  best  indications  that  the  differential  correction  is 
progressing  satisfactorily  is  that  the  range-rate  residuals  from  tracking 
radars  are  small.  Range-rate  residuals  are  most  sensitive  to  timing  errors, 
and  therefore  are  the  last  to  "clean  up”.  Therefore,  a  range-rate  residual 
population  with  -1,  -2  and  -3  exponents  normally  indicates  satisfactory 
convergence. 

For  the  data  environment  available  to  the  analyst  in  the  conduct 
of  decay  exercises,  sophisticated  data  editors  are  not  practical.  Manual 
data  editing,  based  upon  examination  of  the  population  for  "outliers"  is 
often  necessary,  since  the  threshold  editor  used  in  the  program  is  not 
particularly  discriminating.  Avoid  the  temptation  to  remove  whole  passes 
with  self-consistent  large  residuals,  particularly  if  these  fall  in  a  "U" 
region  with  relatively  little  data,  for  the  large  residuals  may  be  sympto¬ 
matic  of  poor  data  distribution  and  editing  may  compound  the  problem. 

4.5  ACQUISITION  MESSAGES 

Look  angles  are  normally  provided  by  the  GLASGP  program  .  During 
the  last  two  to  three  days  of  satellite  lifetime,  however,  these  look  angles 
deteriorate  as  the  SGPDC  theory  fails  to  represent  the  rapidly-changing 
situation  during  decay.  The  recommended  alternative  during  this  period 
is  to  generate  look  angles  with  XYZLA,  utilizing  a  binary  ephemeris  tape 
produced  by  Spiral  Decay.  Experience  has  demonstrated  that  this  procedure 
results  in  look  angles  valid  to  a  few  seconds  in  time,  even  into  the  decay 
revolution. 
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